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We report a measurement of the CP- violation parameter sin 2(3 with B° — > J/ipKg decays in 
which the J/V> decays to hadrons or to muons that do not satisfy our standard identification criteria. 
With a sample of 88 million BB events collected by the BABAR detector at the PEP-II asymmetric- 
energy B factory at SLAC, we reconstruct 100 ± 17 such events, with J/tp — > ir + -K~ir° being the 
most prevalent, and measure sin 2/3 = 1.56 ± 0.42(stat.) ± 0.21(syst.). 
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PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 



Measurement of CP violation in the B-meson system, 
particularly in b — ► ccs transitions, has been a primary 
goal of the BABAR experiment. In the Standard Model, 
these decays exhibit a CP asymmetry that is proportional 
to sin 2/3, where /3 is defined as arg[— V cd Vl/V u V£\, with 



Via the elements of the CKM matrix 



The current 



world average value of sin2/3 is 0.731 ± 0.056 , with the 
B factories (BABAR at SLAC and Belle at KEK) provid- 
ing the most precise measurements 0, Q . The dominant 
decay mode in these measurements is B° — > J/ipK%, 
where only leptonic decays of the J/ip are considered. 
Leptonic decay modes have the advantage of low back- 
grounds, but account for only 12% of J/?/ 1 decays 0- 
Since the current measurements of sin 2/3 are statisti- 
cally limited, in this article we extend the measurement 
through the use of hadronic J/ip decays, as well as pre- 
viously unused muonic decays. 

At the B factories, B° mesons are produced via 
e + e~ — > T(4S) — ► B°B°. For B° mesons produced in 
this manner and decaying to the CP eigenstate J/ipKg, 
sin2/3 appears as the amplitude of a time-dependent CP 
asymmetry. The Standard Model predicts the decay rate 

e -|At|/T B o 

f±(At) = [1 ± sin2/3sm(Am d At)l, 

4r B o 

where the plus (minus) sign indicates that the other, 
"tagging", B° meson in the event decays as a B°(B°), 
At is the decay time of the CP-eigenstate B° meson mi- 
nus the decay time of the tagging B° meson, t b o is the 
B° lifetime, and Ama is the mass difference between the 
two mass-eigenstate neutral B mesons (Am^ is also the 
B° — B° oscillation frequency). The time-dependent CP 
asymmetry is: 



A 



CP 



f+(At) - f-(At) 
/+(At) + /_(At) 



sin2/3 sm(ArridAt). 



Measurement of Aqp requires that a sample of B° mesons 
decaying to J/ipKg be reconstructed, that the flavor of 
the other B° meson in the event be determined, and that 
At be measured. 

A sample of 88 ± 1 million BB events recorded by the 
BABAR detector 5] was used in this analysis. The in- 
nermost component of BABAR is a five-layer double-sided 
silicon microstrip vertex detector with 90° stereo angle, 
allowing precise reconstruction of the location of the B° 
decay vertices along the beam direction. Since the T(4S) 
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is boosted along the beam direction, the difference in po- 
sition between the B° decay vertices in this direction al- 
lows one to measure At. The primary tracking device is 
a 40-layer drift chamber operated with a helium-based 
gas mixture to minimize multiple scattering. The drift 
chamber is surrounded by a Cherenkov particle identifi- 
cation device, and a CsI(Tl) calorimeter. All of the above 
detectors reside in a 1.5 T field generated by a supercon- 
ducting solenoid. The flux is returned via layers of steel 
interleaved with active detectors for the identification of 
muons and detection of neutral hadrons. 

Two types of Monte Carlo simulated events are used 
in the analysis. One, called "full MC" , consists of events 
that are generated according to the known physics of BB 
and continuum production, passed through a detailed 
model of the detector response 0, and reconstructed in 
the same manner as the data sample. The second, called 
"parametrized MC" , consists of events for which the rel- 
evant parameters are randomly generated according to 
the distributions observed in data or in detailed simu- 
lations. For any study where an accurate model of the 
physics or detector response is required, full MC is used. 
Parametrized MC, which can be generated more quickly 
than full MC, is only used to explore the statistical prop- 
erties of the extraction of sin 2/3. 

While many J/ip decays to exclusive hadronic final 
states have been observed 0, the sum of their mea- 
sured branching fractions is less than 20%. To allow 
for the possibility of observing a signal in previously un- 
measured decay modes, we take an inclusive approach in 
the first stage of event selection. Charged tracks are as- 
signed either the electron, muon, pion, kaon, or proton 
mass based on particle identification information, and 
candidates for ir° — > 77 and r\ — > 77 or 7r + 7r~7r° are 
formed. All neutral combinations of up to six tracks and 
neutral mesons are considered (a maximum of two neu- 
tral mesons is allowed) , and those consistent with baryon 
number conservation, strangeness conservation, and Bosc 
symmetry, and having invariant mass mj/^ in the range 
2.80— 3.20GeV/c 2 , are retained for further analysis. De- 
cay modes of the type J/ip — ► KKir are excluded to en- 
sure that the selected sample is independent of the sam- 
ple used in BABAR's previous measurement of sin2/3 0, 
which included B° — ► r) c K^ events with ij c — ► KKir. 

We form K® candidates from a pair of oppositely- 
charged tracks that have invariant mass between 489 and 
507 MeV and a vertex displaced by at least 1 mm from 
the J/ip candidate's vertex. The selected J/ip and K® 
candidates are combined to form B° candidates. Two 
kinematic variables are used to isolate the B meson sig- 
nal: the difference AE between the energy of the re- 
constructed B candidate and the beam energy in the 
center-of-mass frame, and the beam-energy substituted 
mass mES = v^bcam — P*b > where P*b 1S the momentum 
of the reconstructed B and E£ is the beam energy, 
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both in the center-of-mass frame. The small variations 
of E£ eam within the data sample are taken into account 
when calculating ttt-es- Signal events will have AE close 
to and values of ibes close to the B° meson mass. 
Candidates are required to have ibes > 5.20GeV/c 2 and 
\AE\ < 55MeV if the J/ip decays entirely to charged 
particles, and < 105 MeV if the decay includes one or 
more neutral hadrons. The AE selection accepts candi- 
dates within 3cr of the distribution observed in simulated 
signal events. The resolution in toes is 3 MeV, so the 
selection admits a large region at low toes hi addition to 
the region populated by signal candidates. Inclusion of 
this sideband region allows the magnitude of the combi- 
natoric background to be measured. 

Backgrounds arise both from continuum qq production 
and from B meson decays to other modes. The contin- 
uum events tend to have a two-jet topology, in contrast 
to the more spherically symmetric BB events. A set of 18 
variables (described in |7j) that are sensitive to this dif- 
ference are combined in a Fisher discriminant J 7 , which 
is defined to have an average value of 1 for signal and 
-1 for continuum events. The weight of each variable in 
the discriminant is calculated by maximizing the sepa- 
ration between a sample of data taken below the BB 
threshold (and thus composed entirely of continuum qq 
events) and a sample of simulated signal events. We place 
progressively tighter requirements on T as the candidate 
J/ip decay multiplicity increases: for two-body decays 
we require T > —1.14, for three-body decays we require 
T > —0.70, and for higher-multiplicity decays we require 
T > -0.37. 

For three-body J /if) decays, which have a larger com- 
binatoric background than two-body decays, additional 
separation between signal and background is attained by 
considering the angle 9d between the normal to the plane 
in which the momenta of the J/tp daughter particles lie 
and the K® direction in the J/ip rest frame. Conservation 
of angular momentum requires this variable to be dis- 
tributed as cos 2 0d for Jj if; decays to three pseudoscalars 
(the most common type of three-body decays), while it 
is uniformly distributed for BB backgrounds and peaks 
at cos 9d — for continuum qq backgrounds. We require 
candidates to have | cos 0d\ > 0.55. The selection in cos 0d 
and T was chosen to maximize S/ 'v S + B, where S is the 
expected signal and B the expected background. 

There are two classes of BB backgrounds. The first 
consists of candidates formed from a subset of a given 
B meson's decay products, or from a combination of de- 
cay products from the two B mesons in the event. This 
background and the continuum qq background are hence- 
forth referred to as "combinatoric backgrounds" . They 
have a linearly falling distribution in AE, and their dis- 
tribution in Toes may beparametrized by an empirical 
phase-space distribution 8] (henceforth referred to as the 
ARGUS function): 

A(m E s 5 mo , Carg ) oc m E s y/l — (toes/to ) 2 x 

exp(c arg (l - (toes Mo) 2 )), 



where mo is a cutoff mass set to 5.291 GeV (a typical 
center-of-mass beam energy) and c arg is a fitted para- 
meter. 

The second class of BB background consists of B 
mesons that decay to a topology also allowed for J/ipK®, 
but without a J/ip in the intermediate state. These 
"peaking" backgrounds are dominated by B decays that 
have a charmed meson in the intermediate state, so we 
remove any candidates for which afor D* meson within 
2<T of the nominal mass can be formed from the final-state 
hadrons. Since these backgrounds arise from fully recon- 
structed B° mesons, they have the same distribution in 
toes and AE as the signal. 

Since the branching fractions for many of the modes 
that contribute to the peaking backgrounds are not well- 
measured, we must extract the peaking background mag- 
nitude from the data. We do this by performing a two- 
dimensional unbinned maximum likelihood fit to the Toes 
and rnj/jf, distributions. The likelihood function used is: 

L = (n comh A(m E s; m , c arg ) + (n sig + 7ip Cak )G(TO E s)) X 
(Komb + nl eak )C(mj/^;pi,p 2 ) + n sig G(m,j/^)), 

where n com b is the fitted combinatoric background, n® ea:k . 
is the fitted peaking background, n s j g is the fitted signal, 
A is a normalized ARGUS function, G are normalized 
Gaussians, and C is a normalized second-order Cheby- 
shev polynomial with parameters pi. The mean and 
width of G(toes) are fixed to the values observed in high- 
statistics hadronic £?-decay samples, and the mean and 
width of G(toj/^) are fixed to the values observed in 
our J/ip —>■ [i + piT sample for two-body decay modes, 
and to the values observed in full MC events for higher- 
multiplicity modes. The photon-energy resolution in the 
simulated events is degraded to match that observed in 
data. The additional smearing required is 3% of the mea- 
sured photon energy for photons below 100 MeV, and 
decreases with increasing photon energy (no additional 
smearing is required for photons above 1 GeV). 

The J/ip decay modes for which the measured signal 
magnitude is less than its statistical uncertainty are re- 
moved from the analysis. The surviving modes, and their 
contribution to the signal, are listed in Tabled Note that 
no modes including an rj meson are observed, and also 
that no decays with a multiplicity of greater than three 
are visible above background. 

The observation of 28 candidates in the J/ip — ► ir + ir~ 
channel is inconsistent with our expectation of observing 
about one event given the known branching fraction of 
(1.47 ± 0.23) x 10~ 4 2} for this mode. We interpret the 
excess candidates as J/ip — > decays in which both 

muons fail the standard muon selection criteria. Studies 
using simulated events with muon identification efficien- 
cies measured in data confirm that the observed signal 
magnitude is consistent with the J/ip — ► hypoth- 
esis. Since these events do measure sin2/3, and are in- 
dependent of the events used in our previous measure- 
ments [3j, we retain them for this analysis. 



6 



J/ip decay mode 


Signal 


Peaking Bkg 


Comb. Bkg. 


T 

7T 7T 


28 ±8 


84 ± 17 


206 ± 12 


K + K~ 


5±3 


-1 ± 6 


42 ± 5 


PP 


6 ± 3 


1 ± 6 


34 ± 5 


Total h + h~ 


40 ±9 


86 ± 19 


279 ± 13 


After final selection 


28 ±8 


13 ± 3 


15 ± 3 


7T + 7T _ 7T U 


58 ± 17 


104 ± 29 


652 ± 23 


PpTT° 


11 ±6 


9±9 


77 ± 7 


Total h+h-7v u 


69 ± 18 


113 ±30 


716 ± 22 


After final selection 72 ± 13 


19 ±5 


74 ±8 



TABLE I: Observed B° -> J/ipK.% signal and background. 
The combinatoric backgrounds reported are the integral of 
the fitted ARGUS function in the region jties > 5.27GeV/c 2 . 
Except for the rows labelled 'After final selection" , the num- 
bers are measured prior to application of the final selection 



criteria on m 
only. 



and AE. All uncertainties are statistical 



After np 0ak is determined, the following final selection 
criteria are imposed to improve the purity of the sam- 
ple: We recalculate AE with the J/ip candidate con- 
strained to the nominal mass, and define the result as 
AE C . The resolution in AE C is 11 MeV for two-body 
J ftp decay candidates, and 12 MeV for three-body can- 
didates. For two-body J/tp decay candidates we re- 
quire 3.06 < m JH , < 3.12GeV/c 2 and \AE C \ < 33 MeV, 
and for three-body J/ip decay candidates we require 
3.05 < m J/v , < 3.15 GeV/c 2 and \ AE C \ < 35 MeV. 

The efficiency of this selection for peaking backgrounds 
(£pcak) is estimated using full MC events. We define 
£ P cak as the ratio of the area of the fitted Gaussian in 
ibes after the final selection to the area before the fi- 
nal selection. For two-body decay candidates e p0 ak = 
0.15 ± 0.01(stat.) and for three-body decay candidates 
£poak — 0.17 ± 0.02(stat.). An unbinned maximum like- 
lihood fit to the sum of a Gaussian distribution and 
an ARGUS function is performed on the toes distri- 
butions of the surviving candidates. The integral of 
the ARGUS function measures the combinatoric back- 
ground, while the integral of the Gaussian measures the 
sum of the signal and peaking background. Subtracting 



^pcak 



peak "-peak 







from the latter provides an estimate 



of the signal. The toes distributions are shown in Fig.^ 
and the signal and background magnitudes in the final 
sample are reported in Table [I] 

Once the sample of B° — ► J/ipK® candidates has 
been isolated, the extraction of sin2/3 proceeds in the 
same manner as for BABAR's other recent measure- 
ments Q. Information from the final-state particles re- 
coiling against the J/ipK® meson candidate is used to 
determine whether the other B meson in the event was a 
B° or B° at the time of its decay. This is referred to as 
the flavor "tag" . The variables used for tagging include 
the charge of any high-momentum identified electron or 
muon, the charge of any identified kaon, and the charge of 
a slow pion consistent with arising from D* meson decay. 
The efficiency e and mistag rate w are measured using 
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FIG. 1: 771es distributions for candidates for B° — > J /ipK% 
with the J/ip decaying to (a) two and (b) three particles. 
The dotted line represents the fitted combinatoric background 
distribution. The dashed line represents the total background 
distribution, while the solid line represents the signal plus 
background distribution. 



the data as described below, and reported in Ref . 0] ; the 
overall figure of merit for the flavor-tagging performance, 
e(l-2w) 2 , is (28.1 ±0.7)%. 

The extraction of sin2/3 is done using an unbinned 
maximum likelihood fit to the At distribution of the 
candidate events, where the assumed functional form is 
f±(At) convolved with the resolution of the At measure- 
ment, with the mistag probability taken into account. 
The input to the fit consists of both the signal sam- 
ple and a large sample of fully reconstructed B° de- 
cays to DW+tt-, DW+p-, D(*)+cr and J/i(jK * with 
K°* — > iV + 7r — . The B° flavor is known for these modes, 
so this sample constrains a set of parameters describ- 
ing the flavor-tagging performance and vertex resolution, 
The simultaneous fit takes into account any correlations 
between these parameters and the value of sin2/3. The 
result is: 

sin2^ = 1.56 ±0.42 (stat.) 

The At distribution for flavor-tagged signal events is 
shown in Fig. |3 and the CP asymmetry observed be- 
fore correction for backgrounds and mistag probability is 
shown in Fig. [3] In each case a projection of the best-fit 
model is superimposed. 

As a cross-check, the analysis was repeated using a 
sample of B ± — > J/4>K ± events selected in a manner 
analogous to the CP sample, and with the same J/tp 
decay modes considered. This sample yields an apparent 
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At (vs) 

FIG. 2: At values observed in the B° — > J/ij)Ks candidates. 
The plots show the distribution for events in which the recoil- 
ing B meson is tagged as (a) B° and (b) B°. In each plot the 
solid line represents the result of the maximum likelihood fit, 
and the shaded area the contribution of background. 




FIG. 3: At asymmetry observed before correction for back- 
grounds and mistag probability in B° — + J/ipK® candidates, 
with best-fit asymmetry displayed. 



sin2/3 of — 0.13±0.20 (stat.), consistent with the expected 
null result. 

Systematic uncertainties arise from several sources. In 
performing the fit for sin2/3 it is assumed that the back- 
ground has no CP asymmetry. Since some of the back- 
ground is composed of real B° mesons this may not be 
true. Fitting for sin2/3 on a sample composed of candi- 
dates in the mj/^ or AE C sidebands yields 0.18 ± 0.46. 
The signal sample is then refit with the CP asymme- 
try of the peaking background fixed to the ±lcr limits of 
the measured asymmetry, and the observed variation of 
±0.15 in sin2/? is taken as a systematic uncertainty. 

The next most significant systematic uncertainty arises 
from the estimation of the background magnitudes. 
When the sin2/3 fit is performed, the parameter c arg of the 
ARGUS distribution describing the combinatoric back- 
ground is fixed to the central value determined from fit- 
ting the toes distribution. The sin2/3 fit is repeated with 



this value fixed to its ±ler limits, and the observed vari- 
ation in sin2/? of ±0.13 is taken as a systematic uncer- 
tainty. 

The uncertainty on the peaking background arises from 
several sources, the largest of which is the statistical un- 
certainty on n.p 0ak . The next most significant source is 
uncertainty in £ pca k- We estimate the magnitude of this 
uncertainty by observing the variations in e pca k among 
samples of different simulated B° decay modes. In ad- 
dition, one could define e pca k as the efficiency for any 
candidate with toes > 5.27 GeV to pass the final selec- 
tion, rather than defining it as the ratio of fitted Gaussian 
areas. We take the difference between the two definitions 
as a systematic. The estimate of 7ip eak is also subject to 
uncertainty in the distribution of peaking backgrounds 
in to j , which is modeled as a second-order Chebyshev 



polynomial. The variation in 



peak 



when the order is 



changed by ±1 is propagated to the systematic uncer- 
tainty. The accuracy of the fit used to extract the signal 
is verified using background-only samples, such as data 
recorded below the BB threshold or samples of candi- 
dates reconstructed in modes not accessible to the J/ip- 
No statistically-significant signal yields are reported in 
fits to these samples. We assign the largest artificial 
signal yield consistent with these tests as a systematic 
uncertainty. The final source is the uncertainty on the 
resolution of the J /ip peak (which is held fixed in the fit 
that determines Hp 0ak ). Variation of this assumed width 
between values observed in different decay modes yields 
a variation in ".p Cak . The sum in quadrature of all these 
effects totals 25% of the magnitude of n pca k. Repeating 
the fit on many samples of parametrized MC events, each 
of which has the same size and background as the sample 
observed in data, shows that the variation in sin2/3 result- 
ing from a 25% uncertainty in the peaking background is 
±0.07. 

There are potentially differences in the flavor-tagging 
performance and vertex resolution between events with 
hadronic J/?p decays and the other fully-reconstructed 
B decays used to measure these parameters. Perform- 
ing a sin2/3 fit to a large sample of full MC signal events 
with J '/ip — ► 7r + 7T~7r with the flavor tagging and vertex 
resolution fixed to the measured values yields a result 
consistent with the generated value. The statistical un- 
certainty of the result (±0.04) is taken as a systematic 
uncertainty. 

Another systematic uncertainty arises from events in 
which one or more of the final state particles assigned to 
the reconstructed B° in fact originated from the other 
B° in the event. The fraction of such events is negligible 
for two-body J/ip decays, and about 5% for three-body 
decays. Performing sin2/3 fits on full MC samples with 
and without including the incorrectly reconstructed can- 
didates yields a variation of ±0.01 in sin2/3. 

Finally we take into account all the sources of system- 
atic uncertainty that apply to BABAR's previous mea- 
surements of sin2/3 0, except for those specific to the 
B° — > J/ipK® mode, that have not already been specifi- 
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Source 


Uncertainty 


Peaking background CP 


0.15 


Combinatoric background magnitude 


0.13 


Peaking background magnitude 


0.07 


Tagging and vertexing differences 


0.04 


Common to leptonic modes 


0.03 


Misreconstructed signal 


0.01 


Total 


0.21 



TABLE II: Summary of systematic uncertainties on the mea- 
surement of sin2/3. 

cally addressed here. These uncertainties primarily arise 
from limits on our understanding of flavor-tagging and 
vertex reconstruction performance, and yield a variation 
of ±0.03 in sin2/3. 

The systematic uncertainties are summarized in Ta- 
ble [H] The sum in quadrature of all contributions is 
0.21. 

The value of sin2/3 reported in this analysis is higher 
than the world average value of 0.731 ± 0.056. To esti- 
mate the consistency of this result with the world aver- 
age, 10,000 parametrized MC samples with the same sig- 
nal and background magnitudes as observed in the data 
were generated with a true sin2/3 of 0.731. To simulate 
the systematic uncertainty in this analysis and the total 
uncertainty on the world average, a random number with 
Gaussian distribution and a = 0.22 is added to the sin2/3 
result for each sample. Of the 10,000 samples, 629 fluc- 
tuated to a value of 1.56 or greater, indicating that the 
probability of such a fluctuation is 6.3%. 

In summary, we have extended BABAR' s previous sin2/3 
measurement by including J/ipK® modes where the J/tp 



decays to hadronic final states. The result is 

sin2/3 = 1.56 ± 0.42 (stat.) ± 0.21 (syst.). 

Although we searched for many hadronic J/ip decay 
modes, signals were observed only in modes that have 
been previously seen 2J. Further, only in hadron multi- 
plicities of two and three was it possible to observe a sig- 
nal above background. Extending the analysis to the Xc 
and ip{2S) mass regions does not yield additional signifi- 
cant signals, nor is an r/ c signal observed after elimination 
of KKir modes. 
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